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Abstract 

We derive a finite set of nonlinear integral equations (NLIE) for the thermody- 
namics of the one-dimensional s/(4)-symmetric Uimin-Sutherland model. Our NLIE 
can be evaluated numerically for arbitrary finite temperature and chemical poten- 
tials. We recover the NLIE for s/(3) as a limiting case. In comparison to other 
recently derived NLIE, the evaluation at low temperature poses no problem in our 
formulation. The model shows a rich ground-state phase diagram. We obtain the 
critical fields from the T — > limit of our NLIE. As an example for the application 
of the NLIE, we give numerical results for the SU (4) spin-orbital model. The mag- 
netic susceptibility shows divergences at critical fields in the low-temperature limit 
and logarithmic singularities for zero magnetic field. 

PACS: 02.30.1k, 05.70.-a, 75.10.Jm 

1 Introduction 

Since Bethe's seminal solution of the one-dimensional spin- 1/2 Heisenberg chain [H, 
many integrable, natural generalisations of this model have been treated using basi- 
cally the same ansatz. Among them is a multi-component, higher-rank generalisation 
of the Heisenberg chain, first proposed by Uimin for the case of three components 0. 
Later, Sutherland introduced and solved the model for an arbitrary number of parti- 
cle types 0]. The two-dimensional classical model associated with the one-dimensional 
Uimin-Sutherland (US) model is the Perk-Schultz (PS) model [4|]. The case of higher-rank 
representations of the underlying symmetry algebra was treated by Andrei and Johan- 
nesson 0,0]- Affleck calculated the critical behaviour based on non-Abelian bosonisation 
and conformal field theory 0, S] . 
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The traditional thermodynamic Bethe ansatz (TBA) allows for the treatment of 
finite-temperature properties 0, 10, H,[l2]. It uses the string hypothesis and typically 
yields an infinite set of nonlinear integral equations (NLIE). The numerical solution of 
these equations poses a problem as some kind of truncation scheme is necessary. Only in 
the limit T — > one obtains a finite set of equations. Using the TBA approach for the 
general US model [13( , the low- field asymptotics of the susceptibility [LJ] as well as the 



low-temperature asymptotics of the specific heat [15|, [l 



With the help of the quantum transfer matrix ( 



have been derived analytically. 




formalism 17], LL8( and the 



fusion hierarchy of transfer matrices [13 . |2Q . |21| . |22j, 1231 . 12411 . it is possible to rederive 
the TBA equations without using the string hypothesis [25[]. Using the QTM it is also 
possible to derive an alternative set of NLIE consisting only of a finite number of unknown 
auxiliary functions 



26l . 1271 . 128l . 129] , thus allowing for a precise numerical treatment at 



arbitrary finite temperature. Later it was realised that these auxiliary functions provide a 
natural way to exactly truncate the TBA equations [5cJ]. Nevertheless, no straightforward 
way of getting the required auxiliary functions is known. Up to now, NLIE of this type 



31, 32 



34j | . which allows for the generalisation to 



have only been derived for three components at most 

There exists yet another type of NLIE 
an arbitrary number of components [3^]. Unfortunately, these NLIE prove difficult to 
evaluate at low temperature. Instead, high-temperature expansions (HTE) have been 



obtained up to high order [35|, [36(] . Very recently, these equations have been further 



generalised to treat C/ g (sZ(r|s))-symmetric PS models [3 

In this paper, we treat the four-state, s/(4)-symmetric US model in the spirit of [2£ 
We define 14 suitable auxiliary functions, from which we derive a set of well-posed NLIE 
that are valid for arbitrary finite temperature and chemical potentials. The auxiliary 
functions are connected to the fundamental representations of si (4). Our NLIE are the 
natural generalisation of those obtained for the s/(3)-symmetric case (33| . 

The sl(4)-symmetric US model has many interesting applications. As an example, 
we treat the SU(A) spin-orbital model [3^. 139]. The thermodynamic properties of this 
model have already been studied numerically using various methods [37l . 

Eft SHU- The 

ground-state phase diagram in dependence of the magnetic field and the orbital Lande 
factor has also been obtained [43J. In comparison to these methods, we are able to obtain 
highly accurate numerical results for low finite temperatures in the thermodynamic limit. 
Other possible applications of the s/(4)-symmetric US model include an integrable two- 
leg spin ladder system [44j, which has recently been studied numerically using the HTE 
and TBA methods [II, ki]. 

The outline of this paper is as follows. In Section [21 we briefly introduce the (/-state 
US model and show how its thermodynamic properties can be obtained using the QTM 
approach. In Section [3l we concentrate on the s/(4)-symmetric case. We present a set 
of 14 well-posed auxiliary functions, from which we derive a closed set of NLIE. We 
show how the largest eigenvalue of the QTM can be extracted from these auxiliary 
functions. In Section 0J we treat two limiting cases of our NLIE. First, the NLIE and 
auxiliary functions of the s/(3)-symmetric case are recovered by freezing out one of the 
states. Second, the limit T — ► yields linearised integral equations, which are equivalent 
to the corresponding TBA equations and allow for the derivation of the critical fields. 
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R%,{v) = 

a 

Figure 1: Graphical depiction of the i?-matrix as defined in equation 

In Section \5\ we deal with the numerical solution of our NLIE. We briefly introduce 
the SU (4) spin-orbital model as an application and give results for various physical 
quantities. In SectionEl we give a summary of our work and an outlook on open problems. 
Appendix [X] is devoted to details concerning the derivation of the NLIE. 

2 QTM approach to the US model 

In order to fix notation, we begin with a short review of the Uimin-Sutherland (US) 
model Consider a one-dimensional lattice with L sites, where a g-state spin variable 
Oj is assigned to each site j. Each spin a has its own grading e a = (— 1)p(°) = ±1. The 
Hamiltonian of the US model is then given by 

L 
3=1 

where the local interaction operator 7Tj i.fi permutes neighbouring spins on the lattice 
with respect to their grading, 

7r iii+ i|ai . . . ajaj+i ...a L ) = (-If ( a M°«j+i)\ ai . . . a j+1 aj ...a L ), (2) 

and periodic boundary conditions are imposed. The model shows sl(r\s) symmetry, where 
r and s are the total number of states with positive and negative grading (q = r + s), 
respectively. It is therefore a higher rank generalisation of the spin-1/2 Heisenberg chain, 
which is contained as the special case q = 2 and e\ = ei = +1. 

The one-dimensional US model is known to be exactly solvable by Bethe ansatz (B A) . 
The two-dimensional classical counterpart is given by the Perk-Schultz (PS) model [J], 
which is defined on a square lattice with L x N sites and periodic boundary conditions 
in both directions. Variables taking on integer values from 1 to q are assigned to each 
bond of the lattice, and a Boltzmann weight depending on a spectral parameter v, 

<» = <W<W + v ■ (-l) p{a)pM ■ SapS^ , (3) 

is associated to every local vertex configuration (a, /3, //, u), see Figure [TJ We define the 
row-to-row transfer matrix 
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As the .R-matrix (JH]) is a solution to the Yang-Baxter equation, these transfer matrices 
form a commuting family [T(v),T(v')] = for all v,v' E C. Making use of Baxter's 



formula 47]], one recovers the Hamiltonian of the US model from the transfer matrix of 



the PS model at the shift point v = 0, 



d 



H = —lnT(v 

av 



v=0 j=1 



Without breaking integrability, we may add external field terms, 

L q 

H = Hq+ H cy ± = Hq — ^2 ^" n J> ' (6) 

j=l a=l 

where fi a is some general chemical potential associated with state a and the operator 
nj : a counts the number of particles of type a sitting on site j. 

We are interested in the thermodynamics of the US model. Hence, we want to estab- 
lish some connection between its partition function and the transfer matrix of the PS 
model. We therefore consider a second 72-matrix, namely R(v), obtained by rotating the 
graphical depiction of R(v ) clockwise by 90 degrees, 

RZ(v) = Rf a (v) . (7) 

We define the transfer matrix T(v) as the product of matrices R{v) in analogy to (j3|). 
Equation © now also applies to T(v), and as a consequence the relation 

T(-(3/N)T{-f3/N) = e -WNYHo+o(WN)*) ( 8 ) 

is valid for arbitrary inverse temperature (5 and a sufficiently large even integer Trotter 
number N. The partition function of the one-dimensional US model is then given by 



Z = Tve~ m = lim Tr 

N—>oo 



T{u)T{u)) N/2 e - /JW - rt l , (9) 



where u = —(i/N and the traces are taken in the g^-dimensional space. Obviously, the 
partition function of the one-dimensional US model is equal to the partition function of 
a staggered two-dimensional PS model, where the external field can be incorporated by 
modifying the boundary conditions in the Trotter direction. Equation ([9]) is still difficult 
to evaluate, e.g. all eigenstates have to be taken into account. To avoid this problem, it is 
better to consider the column-to-column transfer matrix of the staggered vertex model, 
which is called the quantum transfer matrix (QTM) [l7l. In order to write down the 
QTM in a convenient way, we first define the matrix R(v), which we get by rotating 
R(v) counterclockwise by 90 degrees and changing the sign of the spectral parameter, 

R%(v) = R^(-v) . (10) 
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Then the QTM takes the form 

N/2 

( TQ ™t («) = E e ^ IT R^r^iiv + u)2^ + > - u) , (11) 
M i=i 

where we have introduced a new spectral parameter v, so that the QTMs for all v, v' E C 
form a commuting family, 

[r^™(v),r^™(v')]=o. (12) 

This allows for the diagonalisation by use of the BA. In the end, we are only interested 
in the case v = as the partition function of the one-dimensional US model in terms of 
the QTM is 

Z = Jim Tr(TQ™(0)) L . (13) 

In the thermodynamic limit {L — > oo) one finds that the thermodynamics of the US 
model solely depends on the unique largest eigenvalue of the QTM 48], |49|]. For the free 
energy per unit length, we finally get 

/ = - lim i-lnZ = -ilnA max (0), (14) 
L^oo Lp p 

where A max (v) is the largest eigenvalue of the QTM. 

As noted before, the QTM can be diagonalised via BA. The result for the eigenvalue 
is 

q 

AM = E A »' ( 15 ) 

3=1 

where 

For convenience, we have defined the functions <f>±(v) = (v db iu)^/ 2 and 



(f>-(v) for j = 

n^Li(«-<.) for j = l,...,g-l , (17) 
>+(«) forj = g 



where the complex parameters «•? . are the so-called B A roots and Mj is the total number 
of BA roots in set j. The BA roots have to fulfil the BA equations 

Xj(v{.) 

-1, (18) 



to ensure that all potential poles in the expression (|15|) for A(v), which has to be a 
polynomial of degree N, cancel. The BA equations form a system of coupled nonlinear 
equations for the unknown BA roots. 
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3 Nonlinear integral equations for the sl(4) case 



Let us now turn to the special case of the s/(4)-symmetric Uimin-Sutherland model 
(q = 4, €j = +1 for all j). The largest eigenvalue of the QTM lies in the sector Mj = N/2 



for all j. Of course, it is in principle possible to solve the BA equations (|18p for some 
fixed Trotter number N. But this approach is possible only for finite N and is also quite 
cumbersome to do numerically. As we are interested in the limit N —* oo for deriving 
the free energy of the model, we have to encode the BA equations into a form for which 
this limit can be taken analytically. 

We start by defining some suitable auxiliary functions, which will in the end turn out 
to fulfil certain nonlinear integral equations. For convenience, we wil l use an abbreviated 
notation utilising the Yangian analogue of Young tableaux [ic" 
the function \j(v), we will write a box filled with the letter j, 



5l|, [52, |53|. Instead of 



(19) 



This corresponds to a Young tableau belonging to a vector of the first, four-dimensional 
fundamental representation of s/(4). We also define a Young tableau belonging to the 
second, six-dimensional representation, 



X j {v-i/2)X k (v + i/2), 



(20) 



and for the conjugate four-dimensional representation we have 



\j(y - i)\ k (v)Xi(v + i) 



(21) 



From fusion hierarchy [22j, l25f] one knows that the eigenvalues of the QTMs belonging 
to the three fundamental representations can be written as 



A«(x) 
A^(x) 



[T] + [2] + [3] + [4 

T 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



(22a) 
(22b) 

(22c) 



where the superscripts denote the representations. In all three cases the BA equa- 
tions (fTS|) ensure that the eigenvalues are free of poles. We note that (|22ap is equiv- 
alent to (I15p . We define the following four auxiliary functions for the first fundamental 
representation: 



b[ 1} (x) 



2 + 3 + 4 



(23a) 



v=x+i/2 



6 



&?>(*) 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



v=x— i/2 



We have six auxiliary functions for the second fundamental representation: 

bf\x) = 



+ 



+ 



+ 



+ 



v=x+i/2 



4 2) (*) 



bf\x) 



bf\x) 



+ 



+ 



A 



+ 



+ 



v=x+i/2 



T 



+ 



+ 



+ 



+ 



bf\x) 



bf\x) 



+ 



+ 



+ 



+ 



v=x— i/2 



+ 



+ 



+ 



v=x— i/2 



(23b) 

(23c) 
(23d) 



(24a) 



(24b) 



(24c) 



(24d) 



(24e) 



(24f) 
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And finally, the four auxiliary functions for the third fundamental representation are: 



bf\x) 



4 3) (x) 



bf\x) 





1 






2 






3 




1 




1 




2 


2 


+ 


3 


+ 


3 


4 




4 




4 



(25a) 



v=x+i/2 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



(25b) 



(25c) 





2 






CO 






4 




1 




1 




1 


2 


+ 


2 


+ 


3 


CO 




4 




4 



(25d) 



v=x-i/2 



In addition to the auxiliary functions above, we define a second set of functions, namely 
Bf\x) = bf\x) + 1. These can also be written in a form, where only simple sums of 
Young tableaux appear as factors in the numerators and denominators, and all factors 



appearing in the set of functions 



b { rHx) 



also appear in the set of functions B^ n \x). 



We note that all these factors are partial sums of the Young tableaux appearing in the 
eigenvalues (|22p . The set of functions (|23|) is related to (|25p by a conjugation transfor- 
mation. The set (I24p is self-conjugate in this sense. Furthermore, all auxiliary functions 
are rational functions in terms of the spectral parameter x and are analytic, non-zero 
and have constant asymptotics (ANZC) in a strip —1/2 < Q(x) < 1/2. 

Now, the actual calculation is rather straightforward but lengthy. Therefore, details 
are deferred to Appendix [A] The key idea is to apply a Fourier transform to the loga- 
rithmic derivative of all auxiliary functions and to exploit their analyticity properties in 
Fourier space. We like to stress that, although we are working with arbitrary finite N 
throughout the derivation, in the end the limit TV — ► oo can be taken analytically. 

We finally arrive at a system of coupled nonlinear integral equations (NLIE) of the 
form 

b(x) = —j3e(x) — [K * B] (x) , (26) 



where we have defined 



In b 



(i) 



Mb { l\\nb\ 



(2) 



,ln4 2) ,ln6' 



(3) 



(27) 
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B 



= (inB^, . . .,lnB^\lnB{ 2 \. . . , In B% } , In b[ 3) , . . . ,lnBf } ) 



£ (1) .(2) 



«=(««. 

Convolutions are denoted by 

[/ * g] (x) 



c (2) (3) 



(3) 



f(x-y)g(y) 



dy 
2vr 



(28) 
(29) 

(30) 



The kernel matrix K(x) is a 14 by 14 matrix. As this is too large to be displayed as a 
whole, we divide the matrix into blocks connecting the auxiliary functions from different 
representations, 

(K^\x) K^ 2 \x) KW(x)\ 
K(x)=\K^ 1 \x) K^ 2 \x) K^ 3 \x) . (31) 
\K( 3 ' 1 )(x) K( 3 ' 2 )(x) K( 3 ' 3 )(x)/ 

Obviously, K^'^x) is a 4 by 4 matrix, while K^' 2 ^(a;) is a 4 by 6 matrix, etc. K(x) is 
Hermitian and invariant under reflection along the anti-diagonal, 



K(x) = [K(x)] f , 
Therefore we only need to consider 



K^\x) = 



[K(x)] i>fc = [K(x)] 



15-fc,15-i 



(32) 



K( 1 ' 2 )(a 



K(i.3)( a 



(K Q [x) 


Ki{x) 


K±{x) 


Ki(x)\ 










K 2 (x) 


K (x) 


K^x) 


K^x) 










K 2 {x) 


K 2 (x) 


K (x) 


K^x) 










\K 2 (x) 


K 2 (x) 


K 2 (x) 


K (x)J 










/K u (x) 


K u {x) 


K u (x 


) K 12 (x) 


K 12 (x 


) K 12 (x)\ 


K u (x) 


K u {x) 


K\ 4 {x 


) K ±1 ( 


x) 


K u (x 


) K 12 ( 


x) 


K 13 (x) 


K u (x) 


K\ 4 {x 


) K U ( 


x) 


K 14 (x 


) K n ( 


x) 


\K 13 (x) 


K 13 (x) 


K u {x 


) K 13 (x) 


K u (x 


) K u (x)J 


(K 15 (x) 


K 15 (x) 


Ki 5 (x 


) K 16 (x)\ 








K 15 (x) 


K 15 (x) 


K 18 (x 


) K 15 ( 


x) 








K 15 (x) 


K 19 (x) 


Ki 5 (x 


) K 15 (x) 








\Kn(x) 


K 15 (x) 


Ki 5 (x 


) K 15 (x)J 








(K 3 (x) 


K 4 (x) 


K 4 {x) 


K 4 (x) 


K 4 (x) 


K 6 {x)\ 




K 5 (x) 


K 3 (x) 


K 4 (x) 


K 4 (x) 


K 8 (x) 


K 4 (x) 




K 5 (x) 


K 5 (x) 


K 3 {x) 


K w (x 


) K A (x) 


K 4 {x) 




K 5 (x) 


K 5 (x) 


K w (x) 


K 3 (x) 


K 4 (x) 


K 4 (x) 




K 5 (x) 


K 9 (x) 


K 5 (x) 


K 5 (x) 


K 3 (x) 


K 4 (x) 




\K 7 (x) 


K 5 (x) 


K 5 (x) 


K 5 (x) 


K 5 (x) 


K 3 {x)J 





(33a) 



(33b) 



(33c) 



K( 2 > 2 )( a 



The remaining matrices easily follow from the relations 

K( 3 > 3 )(x) = K^\x) , K^\x) = [K (1 ^(x)]\ 



(33d) 



(34a) 
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K^ 1 ) (x) = [K^ 1 ' 3 ) (x)] 1 , K (3 ' 2) (x) = [K^ (x)] f , (34b) 

[K^) (x)] . fe=[ Ka- 2 )(x)] 7 _ fci5 _.. (34c) 

The kernels are denned as Kj(x) = Kj(k)e lkx dk, where 



k (k) 






k x {k) 


= Kf 4 1 ] ' 1) (fc)+e- fc / 2 -l fc l/ 2 , 


(35a) 


K 2 (k) 


= k^ 1 \k) + ^ 2 -\ k \' 2 , 




km 


= K$ 2) (k), 


(35b) 


k A {k) 


= k {2 f\k)+e- k ' 2 -\ k \l\ 




k 5 (k) 


= kl 2 f ) (k)+e k / 2 ^ 2 , 


(35c) 


K 6 (k) 


= k$ 2 \k) + e- k ~W, 




k 7 {k) 


= 4 2 < 2) (£:)+e fc -l fc l, 


(35d) 


K 8 (k) 


= K {2 f(k) + 2e- k / 2 -\ k \l 2 , 




Kg(k) 


Q\ I /r-i 1 I 1 /r-i 

= K (2 f\k) + 2z k / 2 ~\ k \l 2 , 


(35e) 


K w (k) 


= K {2 f\k) + e~\ k \ , 




k u (k) 


= K { ^\k), 


(35f) 


K l2 {k) 


= k$ 2 \k) + e- k -W 2 -e 


-k/2 

•> 


k 13 (k) 


= k^ 2) (k)+e k -W 2 -e k / 2 , 


(35g) 


k u {k) 


= K$ 2 \k) + e-W 2 , 




k 15 (k) 


= K^f(k), 


(35h) 


k m {k) 


= K [ ( 4 1 ] ' 3) (fe)+e^/^l fc l/ 2 - 




k 17 (k) 


= k^\k)+^ k / 2 -\ k \/ 2 -e k . 


, (35i) 


k 18 (k) 


= k^ 3) (k)+e- k / 2 -^ 2 - 


1, 


k 19 (k) 


= k\^\k)+e k l 2 -\ k \/ 2 , 


(35j) 



with the function 

j>{n,m) (h s_ | fc | /2 sinh(min(n, m)k/2) sinh([g - max(n, m)]fc/2) 
M [) ~ e S mh(k/2)smh(qk/2) ° nm ' W 

We note that in spectral parameter space all kernels can be written in terms of digamma 
and simple rational functions. Nevertheless our notation is more useful here as the nu- 
merical treatment of the NLIE can be conveniently done in Fourier space. The functions 
Kfa\' m \ x ) = f^oo K^' m \k)e lkx dk are related to the S'-matrix of elementary excita- 
tions [54| via 



The bare energies in (f29|) are (x) = vffl (x) + Cj , where 

v ffr x ) = 2 JL gimWg) ; (38) 

[ g ] q cosh(27rx/g) — cos(im/q) 

and the constants are given by 

4 = (-3/ii + fi 2 + M3 + M4)/4, 4 = ~ 3 ^2 + M3 + M4)/4, (39a) 
4 1) = (Ml + M2 - 3/i 3 + /i 4 )/4, = (Ml + + M3 - 3/U 4 )/4, (39b) 

(2) (2) 

c\ J = (-Ml - M2 + M3 + M4J/2 , = (-Mi + M2 - M3 + M4J/2 , (39c) 
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(2) (2) 

dj = {-pi + M2 + M3 - M4J/2, q = (Ml - M2 - M3 + M4J/2 , (39d) 

(2) (2) 

4 = (Ml - M2 + M3 - M4j/2 , q = (mi + M2 - M3 - M4)/2 , (39e) 

(3) (3) 

c\ = (-Mi - M2 - M3 + 3m4)/4, 4 = (~Mi - M2 + 3M3 - M4)/4, (39f) 

c 3 = (-Mi + 3M2 - M3 - M4)/4, q ; = (3mi - M2 - M3 - M4)/4. (39g) 

Finally, the largest eigenvalue of the QTM can be written in terms of the auxiliary 
functions, 

1 - 7 - 2 - z E ^ + E E * ln B l n) ] (°) > ( 4 °) 

J=l ' 71=1 j = l 

where d n = ( ) is the dimension of the nth fundamental representation. Therefore, the 
problem of solving the infinitely many BA equations (]18p in the limit N — > 00 has been 
reduced to finding a finite set of functions satisfying the NLIE f|26|) — (|39[) . The NLIE is 
valid for arbitrary finite temperature and chemical potentials. 



4 Analytical investigation 
4.1 Investigation of the sl{3) limit 

First we want to show how our formulation (I26p - (l40p reduces to the known NLIE for the 
sZ(3)-symmetric case [33] by freezing out one of the states. We choose the state a = 4 
and accordingly treat the limit M4 — * —00. We observe that 





;*) 


= 00), 


&«(*) 


= 




b^{x) 


= 




(41a) 




[x) 


= 0(e^ 4 ) , 


bf\x) 


= 




b { 2 2) {x) 


= 




(41b) 


bf 


[x) 


= 0(e^ 4 ) , 


bf{x) 


= 




bf{x) 


= 




(41c) 


bf 


[x) 


= 0(e^ 4 ) , 


bf\x) 


= 




bf{x) 


= 




(41d) 


bf 


[x) 


= 0(e^ 4 ) , 


bf\x) 


= 


e^ 4 ). 








(41e) 



Therefore only seven of the auxiliary functions survive. We can regard 



b { l\x) 



bf\x) 



bf\x) 



b£\x) 



bf{x) 



l(3), 



■A (X) = b<?\x) = 0. 



We also conclude that b^\x)/B\ 

1.(3)/ 



l(3)/ 



(42) 



1. Using this information, the equation for 



In 6^ (x) linearises and can be solved analytically. We get 

Mi + M2 + M3 



liiB[ 3) {x) 



(3 W{x) 



+ M4 



{x) 



v^Khb^b^b^) 



{x), (43) 
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where 



W(x) 



-|fc|/2 



+ l + e A 



i ikx dk . 



(44) 



(2) (2) 

Substituting this into our NLIE and relabelling b 4 (x) to 63 (x), we are again left with 
a NLIE of type (|26p . but with only six remaining auxiliary functions belonging to the 
two fundamental representations of sl(3). Here we get the kernel matrix 



K(i,i) (x ) k(i,2) (x 
K(i,2) (x) k(i,i) (x ; 



K(x) = 

where the submatrices are given by 

K^(x) = 

The Fourier transformed kernels are found to be 



(K Q {x) 


Kt 






K2{x) 


K 






W*) 


K 2 




^o(x) 


(K 3 {x) 


K 3 




K A (x) 


Kz{x) 


K 6 




K 3 (x) 


\K 5 (x) 


K 3 




K 3 (x) 


found to be 







k (k) 






K 2 (k) 


= k^\k) + 


e*/2-|fc|/2 


k A (k) 


= k§ 2 \ k ) + 


e -*-|fc[/2 


k % (k) 




e H*l/2. 



-fc/2 



For the bare energies, we get e^(x) 



Vj^ (x) + with the constants 



(-2^1 + ^2 + ^3)73, 
(Ml + - 2/i 3 )/3, 
-2 - (-Ml + 2^2 -M3)/3, 
The largest eigenvalue is finally given by 



-3 

,(2) 



.(1) 
"2 

.(2) 
-1 

,(2) 



(/il - 2^2 +M3)/3, 

(-Ml - M2 + 2^ 3 )/3- 
(2/xi - /i 2 - M3)/3. 



lnA max (0) = -/? 1 



7T 

3^/3 



In 3 



i=i 



n=l j=l 



(45) 



(46a) 
(46b) 





= k^ 1 \k)+e- k / 2 -^ 2 , 


(47a) 


#3 (k) 




(47b) 


K 5 (k) 


= K^ 2) (k)+e k -^ 2 -e k / 2 , 


(47c) 
(47d) 



(48a) 
(48b) 
(48c) 



(0) . (49) 



As expected, this is exactly the known NLIE for the s/(3)-symmetric case [3c 

We can also recover the explicit form of all auxiliary functions for the s/(3) case. 
We drop all Young tableaux that contain 
X A (v) — > in the limit [i A - 



4 in the auxiliary functions (|23l )— (|25[) as 



-00. We have already seen that seven of the functions 
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(2) (2) 

become zero, while one function diverges. After relabelling b\ (x) to 63 (x) again, the 
remaining six functions take the form 



4%) 



2 + 3 



® 



v=x+i/2 



v=x—i/2 



[2] • as + [2] + in) 



b£\x) 



bf\x) 



, bf\x) 



+ 



+ 



+ 



, (50a) 



(50b) 



v=x+i/2 



+ 



(50c) 



v=x— i/2 



Although only two sets of BA roots are present in the sl(3) case, these auxiliary functions 
still contain the function qs(v), which enters through \s(v). But, as already indicated by 
the definition of (73 (v) for the case q = 3 in (fT7|) . we have to demand (v) = (p+(v) in the 
limit ^4 — » —00 in order that the previously derived NLIE can also be derived directly 
from the auxiliary functions (|50p . The auxiliary functions are equal to those presented 
in @- 

We finally note that our choice of freezing out the state a = 4 is completely arbitrary. 
Choosing one of the other states yields, after relabelling some of the indices, the same 
NLIE and auxiliary functions. For fii — > —00, we find that qi(v) = (p-(v), while q2(v) 
and qs(v) contain the remaining BA roots. For /U a — > —00 with a = 2 or 3, we find that 
q a -i(v) = q a (v). 



4.2 Limit T — > and critical fields 

We divide the NLIE (|26p by (3 and define rescaled auxiliary functions by 



in). 



In the limit T (0 



E { r\x 



-> 00) we get 



Ef\x) 



1 



(x). 



,(») 



(*) 



(x) if 5R(eJ n) (x)) > 
if 5R(eJ n) (x)) < 



(51) 



(52) 



Obviously, auxiliary functions with negative real parts for all x G R do no longer con- 
tribute to the ground-state energy, because e + ^\x) = for these functions. 

Without loss of generality, we choose the chemical potentials to be ordered, fi\ > 
fJ-2 > ^3 > ^4- Changing the order just amounts to some permutation of indices in the 
following calculations. In our case, we observe that the only remaining auxiliary functions 
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are ef~\x) for n = 1,2,3, i.e. one from every representation. The NLIE linearise and 
take the form 

ef\x) = -V$[x) - - £ [K^ m) * e+i m) ] (*) (n = 1, 2, 3) . (53) 

m=l 

It follows that the remaining auxiliary functions are real and symmetric with respect to 
the spectral parameter. The ground-state energy is given by 

fo = 1 " \ ~ f m 2 - \ £ H ~ t * e+ i n) ] (0) • ( 54 ) 

j=l n=l 

These equations have a particularly simple solution if all chemical potentials are equal, 
Ml = = A*3 = M4- In this case, we get e{\x) = —Vffi(x) and e +( [ l \x) = for all n. 
Therefore, the ground-state energy is just / = 1 — 7r/4 — 31n(2)/2. 

In general, depending on certain differences of the chemical potentials, the ground 
state can be in one of four possible phases. We start with the phase, where all degrees 
of freedom are frozen out, i.e. only the state a = 1 survives. In this case, we have 
e\ (x) = e + ^f l \x) for all n. As a consequence, (|53|) can be solved analytically. We find 
the restriction [i\ — fi2 > 4 and obtain 

e W ( x ) = /i! - pt 2 - ^ i ' e i 2) ( x ) = ^2 - ^3 , ef ) (x) = /i 3 - /x 4 , (55) 

while the ground-state energy turns out to be /o = 1 — fix. As expected, the ground state 
is fully polarised. We call the point fj,\ — fj,2 = 4 the first critical field. 

Below the first critical field, i.e. if [ix — fJv < 4, the function e±(x) possesses two 
symmetrically distributed real roots, and we have ef~\x) = e + ^ l \x) only for n = 2,3. 
We can still solve (|53l) for the latter two functions and get 



«£>(*) = -V«(x) + ^ - [*g' l > * e+?>] (x) , (56a) 

eS 2) (x) = K^\x) + i X> - M3 ) - * e + S 1} ] (s) , (56b) 

j'=i 

e[ 3) (x) = fi 3 - (56c) 
For the ground-state energy, we arrive at 

/o = l-21n2-^±^-[^*e^](x). (57) 

This is exactly the T = behaviour of the spin- 1/2 Heisenberg chain. Two states, a = 1 
and 2, are present in the ground state. Note that these equations are valid only above 
the second critical field, i.e. as long as (x) > for all i£i From (l56bl) . we find the 
restriction 

2 

5>i - /i 3 ) > 41n2 + 2 • [Vg } * e + i 1} ] (0) , (58) 
i=i 
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where the positive convolution term unfortunately still depends on the function (x), 
which is not explicitly known. The convolution term vanishes if fi\ = fi2- 

Below the second critical field, the state a = 3 also contributes to the ground state. 
Both e^'(x) and e±\x) possess two real roots, and only (x) = e +( f\x) remains valid. 
Here, we recover the T = behaviour of the s/(3)-symmetric US model, 

4 B >(s) = -V${x) - - [K^ m) * e + i m) ] (x) (n = 1, 2) , (59a) 

m=l 

e?(x) = -W(x) + \ i V lt n) * e+ ^] ^ ' ^ 59b ) 

j=l n=l 

where cf^ are the constants of the sZ(3) case defined in (|48jl . The ground-state energy 
can be calculated by use of 

/o = 1 " A " ^ 3 - \ H ~ E * e + i n) ] (0) . (60) 

j=l n =l 

These equations hold as long as we are above the third and last critical field. We see 
from (I59bh that the restriction is 

3 2 

5>j - M4) > 7rV3 - 3 In 3 + 3 • E [V[|~ n) * e + i n) ] (0) , (61) 

j=l n=l 

where we again have no explicit expression for the positive convolution terms, which 
vanish if [i\ = Hi = 

Finally, below the third and last critical field, all auxiliary functions possess two 
symmetrically distributed real roots. Therefore, equations (f53|) and (f54"l) have to be used 
without further simplification. In this case, all four states contribute to the ground state. 
We note that the equations (|53p and (I54p are equal to those one can get either from the 
traditional TBA equations in the limit T — > or directly from the BA equations for the 
Hamiltonian 13J] . 



5 Numerical investigation 

The NLIE of the s/(4)-symmetric Uimin-Sutherland model (|26p are of a type that allows 
for an efficient numerical solution by iteration and use of the fast Fourier transform (FFT) 
to calculate the convolutions. As initial values, some discretised functions In Bj n \x) 
are taken. The FFT is applied and the right-hand sides of the NLIE are computed in 
Fourier space. After that, the functions In 6^ (x) are obtained using the inverse FFT to 

eventually yield a new approximation for the functions lnBj n \x). These steps have to 

be repeated, each time starting with the previous approximation for lnBj n \x), until the 
numerical error is small enough. Finally, the free energy is calculated using (I14p and (|40p . 
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Prom our NLIE it is also possible to directly calculate derivatives of the free energy 
with respect to some parameter p. We just have to consider the corresponding derivatives 
of the NLIE and of the expression for the eigenvalue. For the first derivative, we exploit 
the relation 

* ^ (I) = ^'* ! W ( } 

and use an additional relation for the calculation of the second derivative, 



5?' nBr (I) = ^ 1^ l*' ni; W J + 5? ln6 > W j • < 63 > 

In this case, the numerical calculation is done step by step. After solving the unmodified 
NLIE to obtain \nBj n \x), its first and finally its second derivative are solved in an 
analogous way. In each step, results from the previous calculations are used. 

In the following, we will provide some numerical results for a special application of the 
one-dimensional s/(4)-symmetric Uimin-Sutherland model. We consider the Hamiltonian 

L L 

H = ^SjSj+i + l/2)(2r i r i+1 + 1/2) - £ (g s hS* + g T hrf) , (64) 
3=1 3=1 

which describes the 5'C/(4)-symmetric case of a SU(2) x SU(2) spin-orbital model at the 
supersymmetric point. Sj is a SU(2) spin-1/2 operator acting on the spins and Tj is a 
SU (2) spin-1/2 operator acting on the orbital pseudo-spin degrees of freedom. We have 
allowed for an external magnetic field h, which couples to the spins and orbital pseudo- 
spins with Lande factors gs and g T , respectively. Clearly, the Hamiltonian is equivalent 
to the sZ(4)-symmetric US Hamiltonian if we use the basis 

|l) = |Mr), |2> = |Mr>, |3) = |Mr), |4) = ||5lr) (65) 

and accordingly set 

Mi = (gs + 9r)h/2, [i 2 = {gs - 9r)h/2 , (66a) 

= ~(gs - 9r)h/2 , (14 = -(g s +g T )h/2. (66b) 

We are mainly interested in the entropy S, specific heat C, magnetisation M and 
magnetic susceptibility x, which are defined by 

df d 2 f df d 2 f 



In Figure 2 on the next page and Figure 3 on page 18\ results are shown for the case 



gs = 1 and g T = 0, for which the magnetic field couples only to the spins. In this case, 
we have (X\ = (i% = h/2 and (13 = (14 = —h/2. Therefore we know from our analytical 
investigation that there is only one critical field exactly at h c = 2 In 2 ~ 1.39. Below 
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T 

(b) Specific heat vs. temperature. 



Figure 2: Entropy and specific heat of the spin-orbital model at gs = 1, g T = for 
various magnetic fields. The insets show the low-temperature parts. The critical field is 
h c = 2 In 2 « 1.39. 
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(a) Magnetisation vs. temperature. 
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(b) Magnetic susceptibility vs. temperature. 

Figure 3: Magnetisation and magnetic susceptibility of the spin-orbital model at gs = 1, 
g T = for various magnetic fields. The critical field is h c = 2 In 2 « 1.39. 
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0.25 pr 




Figure 4: Magnetic susceptibility of the spin-orbital model at h = for temperatures 
down to T = 1CT 10 . The cross denotes the ground-state value x(0) = 2/ir 2 ~ 0.2026. 
The inset shows the low-temperature part of the susceptibility using a logarithmic scale. 



the critical field, all four states contribute to the ground state. Above, the spins are 
fully polarised and only the orbital degrees of freedom remain. Hence, the ground state 
resembles that of the spin-1/2 Heisenberg chain. This phase transition is clearly exposed 
by the numerical data. The low-temperature slopes both of the entropy and the specific 
heat increase from 2 at h = to infinity at h = h c , whereas they are 1/3 for h > h c . 
Moreover, the magnetisation data show the expected saturation behaviour for h > h c . 
Notice also, that the magnetic susceptibility diverges at the critical field. Below, the 
value at T = stays finite; above, it drops to zero. 

The magnetic susceptibility at h = is particularly interesting as it is expected to 
show a characteristic singular behaviour at T = due to logarithmic corrections. Indeed, 
this is confirmed by our results for the low-temperature susceptibility, see Figured! Even 
for the lowest plotted temperature, T = 10 -10 , the susceptibility is still well above the 
ground-state value x(0) = 2/ir 2 . For the spin-1/2 Heisenberg model, these corrections 



57|, [5aj; similar results for sZ(c/)-symmetric 



have already been treated in detail [55|, [5f 
US models are known 

We also consider the spin orbital-model at gs = 1 and g T = 2, which corresponds 
to the spin-3/2 interpretation. Here we have all three possible types of phase transi- 
tions. Numerical results for this case showing the rich resulting structure are plotted 
in Figure 5 on the following page and Figure 6 on page 21 Again, we observe that the 
low-temperature susceptibility shows characteristic singular behaviour at h = and di- 
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(a) Entropy vs. temperature. 
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(b) Specific heat vs. temperature. 

Figure 5: Entropy and specific heat of the spin-orbital model at gs = 1, g T = 2 for 
various magnetic fields. The critical fields are h c \ = 4, h C 2 ~ 0.941 and h c ^ ~ 0.370. 
Note the numbering of the critical fields as discussed on page [22j 
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(a) Magnetisation vs. temperature. 
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(b) Magnetic susceptibility vs. temperature. 



Figure 6: Magnetisation and magnetic susceptibility of the spin-orbital model at gs = 1, 
g T = 2 for various magnetic fields. The critical fields are h c \ = 4, h C 2 ~ 0.941 and 
h C 3 w 0.370. Note the numbering of the critical fields as discussed on the next page. 
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Figure 7: Phase diagram of the spin-orbital model for magnetic field h and Lande factor 
g T , where gs = 1 is held constant. There exist five different phases (I-V), for details see 
the text on the current page. 



verges at the critical fields. The highest critical field is at h c \ = 4, the other two have 
to be calculated numerically as only the lower bounds h C 2 > 41n(2)/3 ~ 0.924 and 
h c3 > vr/(2V3) - ln(3)/2 » 0.358 are known explicitly. We find the remaining critical 
fields to be h C 2 ~ 0.941 and h c3 ~ 0.370. Note the numbering of the critical fields, where 
at each field h C j the number of involved degrees of freedom changes from j to j + 1. 
The advantage of this scheme is that the critical fields and h C 2 also appear in the 
spin-1 interpretation of the s/(3)-symmetric US model 33J], while only h c \ remains in 
the spin-1/2 Heisenberg model. 

Figure [7] shows the complete phase diagram of the spin-orbital model depending on 
the magnetic field h and the orbital Lande factor g T , while the Lande factor of the 
spins is set to gs = 1. Obviously, there exist five different phases. Above h c \ (I), all 
spins and orbitals are fully polarised. Between h c \ and h C 2, we have to distinguish the 
regions g T > 1 (II) and g T < 1 (III). In the first region, the orbitals are fully polarised, 
while the spins are only partially aligned. In the latter case, it is the other way around. 
For g T = 1, we have a direct transition from phase I to phase IV, because h c \ = h C 2- 
For a magnetic field below h C 2, but above h c3 (IV), the spins and orbitals are both 
partially polarised, while the state |!sIt) is still completely suppressed. For h < h c3 (V), 
all possible spin configurations contribute to the ground state. Notice that h c \ tends to 
infinity for g T — ► 0. As we have seen before, only one phase transition survives for g T = 0, 
where the magnetic field couples only to the spins. The phase diagram presented here is 
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qualitatively in perfect agreement with [431 ] . where a finite system of 200 sites has been 
used for the calculation. 

6 Summary and outlook 

We presented a set of suitable auxiliary functions that allowed us to derive a finite set of 
NLIE for the thermodynamics of the s/(4)-symmetric one-dimensional US model. These 
NLIE are well posed for an efficient numerical treatment using the fast Fourier transform 
as they are of convolution type. Moreover, they are valid for the complete temperature 
range and arbitrary chemical potentials. We have provided analytical results for several 
limiting cases of our NLIE. Thus, we recovered the previously known NLIE of the sl(3) 
case and the linearised integral equations for T = 0. Both provide further support for 
the validity of the NLIE. Using the latter limit, we were able to derive all critical fields. 
We also gave some numerical results for the SU (4) spin-orbital model as an example for 
the application of the s/(4)-symmetric US model, achieving high accuracy even at low 
temperatures. 

This is the first time that NLIE of this type were derived for the sZ ((^-symmetric 
US model with q > 3. Unfortunately, up to now there is no general scheme to construct 
the needed auxiliary functions for the case q > 5. The functions presented here were 
basically found by trial and error. Even at the level of the NLIE, it proves difficult to 
generalise the structure of the somewhat complicated kernel matrix. Once the NLIE are 
known for some number of states q, there is a straightforward way to extract the NLIE 
of all sZ (ra)-symmetric cases with n < q by freezing out one or more of the states. Going 
the opposite way fixes some of the structure, but does not provide enough information to 
deduce the entries of the kernel matrix lying on the anti-diagonal. Nevertheless, we hope 
to overcome these difficulties and to generalise the approach to cover the whole class 
of sZ (c/)-symmetric models in the near future. For that, we expect that the number of 
required auxiliary functions will be equal to the sum of the dimensions of all fundamental 
representations of sl(q), i.e. 2 q — 2. We note that a further generalisation of this type 
of NLIE from the fundamental representation to higher rank representations [5[ should 
be quite straightforward. To achieve this, one has to consider slightly modified auxiliary 
functions, which can be used to truncate the TBA equations at an arbitrary step, in 



analogy to |30j]. 



Another open question is the connection between the kernel matrix and the complete 
S'-matrix of elementary excitations [o^ ]. In the T = case, all kernel functions can 



be obtained from corresponding S'-matrix entries 54|, the generalisation to T > is 
unknown. 

The treatment of graded models is a further direction of g eneralisation. Up to now, 



only the sZ(2|l)-symmetric case has been treated this way 3ll. l32l|. Although supersym- 
metric models are generally more difficult, it turns out that, at least for the sZ(2|l) case, 
the corresponding NLIE are even simpler than the NLIE of the si (3) case as only three 
auxiliary functions are needed. The NLIE for the graded q = 4 models, e.g. the s/(2|2) 
case [60J, also seem to be simpler than the NLIE presented here. We hope to report soon 
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on details concerning this model. 
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The auxiliary functions, as defined in f|23|) — (|25 [) . are rational functions. Moreover, we 
have numerically checked that they are analytic, non-zero and have constant asymp- 
totics (ANZC) in a strip —1/2 < < 1/2, which includes the real axis. Therefore, 

we are allowed to apply a Fourier transform to the logarithmic derivative of all auxiliary 
functions, 



In the cases k < and k > 0, we can close the integration path to a contour above 
and below the real axis, respectively. Hence, it is important to know in which of these 
regions the first-order poles of the logarithmic derivatives are located and how they can 
be further classified. 

Let us first introduce some notation. We define the constants 



A Derivation of the NLIE 




(68) 



a 



(l) 



e A"l _|_ e / 3 M2 _|_ e /3M3 _|_ ^ 
e /30l+M2) _|_ g/3(/W+M3) _|_ e /3(m+M4) 



(69a) 



a 



_|_ g/3(/^2+A13) _|_ e /3(/i2+M4) _|_ g/3(M3+M4) ^ 
e /3(W+M2+M3) _j_ g/3(^l+M2+M4) _|_ e /3(A*l+^3+^4) _|_ e (8(M2+A'3+M4) 

e^ 1 + e^ 2 + e^ 3 , x? } = e' 3 '" 2 + e^ 3 + e^ 4 , 

g^(W+M2) _|_ g/3(Ml+M3) _|_ g/3(Ml+M4) _j_ g^(M2+M3) _|_ g/3(^2+M4) ^ 

e /3(/il+A*3) _|_ e /3(/il+M4) _|_ e /3(M2+M3) _|_ e ^(M2+M4) _|_ e /3(^3+M4) ^ 

e /3(/il+M2) _|_ g/3(^l+M3) _|_ e /3(M2+M3) ^ 

g/3(M2+A13) _|_ e /3(/i2+A*4) _|_ g/?(M3+M4) _ 



(69b) 



a 



(3) 



(69c) 
(69d) 
(69e) 
(69f) 
(69g) 
(69h) 
(69i) 




.(1) 



(3) 



(3) 



For the eigenvalues (|22|) . we find the factorisation 



A( 2 ) (x) = 0_ (x + i/2)<£+ (x - i/2) A( 2 ) (x) • a^ 2 ) , 
tS*\x) = 4>^(v)(j)-{v + i)<t>+(v - i)(t) + (v)A (3 Hx) ■ a (3) 



(70a) 
(70b) 
(70c) 
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where all M n \x) are polynomials of degree N with the highest coefficients being one. 
For the other terms in the auxiliary functions, which are generated by sums of certain 
Young tableaux, we also find that several of their potential poles vanish due to the BA 
equations (fTBj) , We can therefore write 



B + L2 



2 + 3 



3 + 4 



\T} + [2] + [3 



2 + 3 + 4 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



92(a) 
cj)^{x)<i) + {x)qf ) (x) 
qi(x)q 3 (x) 

(j}-{x)qf' \x) _ 

92(a) 
<p + {x)X^\x) 

93(a) 
4>-{x)x£\x) 

91(a) 

0-(a+ ^)4>+ (x 



H>2 



X 



(2), 



(2) 



92 0r 



<j>^{x + ±)<t> + {x-$X%\ X ) 



q 2 (x+ £ 

<P-{X+ Tt)<t>+{X- 



Xi 



(2) 

x 2 



.(x + i)Jf®(x) 



(3) 



93(s + k) 



4>-{x+ h)4>+{x 



(3), 



X 



91 (a - 5) 



•xi 



(71a) 
(71b) 
(71c) 
(71d) 
(71e) 
(71f) 

(71g) 
(71h) 
(71i) 



All functions q^\x) and Xj l> (x) are polynomials of degree N, with the exception of 



('2') 1 

X| (x) and A^ (x) being of degree 3N/2, and with a highest coefficient of one. We note 
that the roots of the functions q^ (x) provide additional solutions to the corresponding 
BA equations (|18j) and are called the hole- type solutions. 

The roots of the polynomials can be obtained from numerical solutions of the BA 
equations at finite N . We find that, in the complex plane, all roots are located on slightly 
curved lines close to horizontal axes at certain heights. Each curve contains N/2 many 
roots. The corresponding height values for each polynomial are given in the following 
table: 



.(2) 



qj(x): 
qf\x): ±1 
A{ 2) (x): +1/2, ±3/2 
Xj 3) (x): ±3/2 



A^(x) 
Xf\x) 
Xf\x) 



Hi + 1)/2 
±1 

-1/2, ±3/2 



(72a) 
(72b) 
(72c) 
(72d) 
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We like to stress that the deviations from these axes remain small even for large N. 

Next, we write the auxiliary functions b^\x) in a factorised form, from which one 
can easily read off the locations of all roots and poles, 



x 



x 



X 



X 



X 



X 



X 



X 



X 



X 



X 



X 



X 



X 



4>- 



4>- 



<t>- 



x - h)<t>+{x + i)qi(x + fi) e^ 1 



(f) _( x + l)xi 1 \x+^ 
x - i)<M* + » ) ft (x + §i)Af >(*) e^ X f 



qi (x+^ q f\x + ^){x) 



ip 3 a 



(2) 



h)4>+( x + 3)92(s - §i)<? 3 (x + |i) e 



,2/3//3 



( 2 ), 



93(x-^ 
x- + i)g 3 (x- §i 



(2) 
X 2 



x -\)4> + (x + \)q 2 (x + 2i) e 



— w ' 

Xi 

,/3(Mi+M2) 



4 2) (x+i) 



(2) 
X 2 



2 V- ' 2 

x - i)0+(x + i)qi(x + i)g2(x - i)<? 3 (x + 2i) e 



/3( M i+2 M3 ) 



(/2 (x + i)xf ) (x)A^(z + ! 
x - i)0 + (x + i)gi(x + i)g 3 (x - i) e- 



(3), 



(i) (3) 
Xi X 2 



^(/il+jU4) 



4' l) (x)A( 1 )(x) 



(i) 



x - i)(j>+(x + i)qi(x - 2i)q 3 (x + 2i) e : 



2/3(/i2+M3) 



9 2 



x)A( 3 )(x) 



ip 2 a 



(3) 



i)0 + (x + i)qi(x - 2i)q 2 (x + i)g 3 (x 



q 2 (x-i)xi 1 \x)Xr>(x-^ 



x - i)<j) + (x + i)q 2 (x 



-(3), 
2i) 



3 /3(2^ 2 +^4) 

(1) (3) 
X 2 Xl 



X{ 2) (x 



3 /3(M3+M4) 

(2) 

Xi 



.(s + f i)<? 3 (x + f i) ePtei+w+w) 



(3) 



h)(p + (x + h)q 2 (x + h)q 3 (x - ±) e 



3 /3(/il+/i2+/i4) 



g 3 (x + |i)x} 2) (x) 



(2) 



(x + |i)g a (x-|i)X«( a : + |) e 



qi {x-%i)qi h \x+i)A( 2 )(x) 
%i)<i>+(x + §i)gi(z - fi) e^+^+w) 



c>/3(/il+/i3+i"4)^( 1 ) 
3^2^30(2) 



(73a) 
(73b) 
(73c) 
(73d) 
(73e) 
(73f) 

(73g) 
(73h) 

(73i) 
(73j) 
(73k) 
(731) 
(73m) 
(73n) 
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The uppercase auxiliary functions B^(x) factorise in an analogous way, 



(i) 



B. 



(i) 



B. 



(i) 



B 



(i) 



B 



(2) 



B. 



(2) 



B 



(2) 



B 



(2) 



B 



(2) 



(2) 



(3) 



13. 



(3) 



(3) 



5 



(3) 



0_(x + l)X«(x + i)"x«' 

^(x + i)xf)(x) 
^(x + ^^x + ^At 2 )^) ^ 3fl(2) 

(x + i)X«(x-^) ^ 



(2), 



X 



q 3 (x - |)X 

g 3(3; _ i)A(D(x 

g 2 (x + i)A( 2 )(x + i 



(2) 



,(1) 



(2) 
X 2 



qi h \x)X?>(x+±) 



( 2 ), 



(2) 
<P2X 2 



q 2 (x + i)x[ 1) (x)X^(x + ^ X W X ® 



qf\x)M l ){x) 
X[ 3 \x 



(i) 



(x + i) 



(3) (3) 
Xl X 2 



gf ) (x)A( 3 )(x) 
qf\x)xf\x 



tp 2 a 



q 2 (x-i)xi 1 \x)Xf'(x-^ 
q 2 (x-i)AW(x-± 



(3), 



(3) ' 

(2) 
<P2Xl 

x? } xS 3) 



X 



(2), 



(2) 



g 3 (x + |i)A( 3 )(x + |) 
+ (x + |i)x{ 3) (x) 



,(3) 



(3) 



<B,(x + fi)x[ 



(3), 



,(3) 



(2), 



xf } (x)X 



(3), 



(2) (3) 
Xl X2 



gi (a; 
g-i (a; 



i)C ; (x + |)A(2)(x) ePw V3 a 
i)A] 3) (*-i) a® 



/,(2) 



U)X^(x) 



e ^ lx (3) 



(74a) 
(74b) 
(74c) 
(74d) 
(74e) 
(74f) 
(74g) 
(74h) 
(74i) 
(74j) 
(74k) 
(741) 
(74m) 
(74n) 



Now, we apply (I68|) to all auxiliary functions. For brevity, we just treat the case k < 
here as the calculation is completely analogous for k > 0. In this case, we just have to 
deal with roots and poles, which are located above the real axis. For the functions of 1 (x) 
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we get the result 

feW (k) = e k/2 ^ (k) - e~ k/2 X^ ] (k) , (75a 

6« (k) = e fe / 2 ^_ (k) + X {3) (k) - e- fe / 2 # (k) - A^ (k) , (75b; 

&W(fc) = e fe / 2 0_(A;) + e 3k/2 q 2 (k) - e k /%(k) - xf\k) , (75c 

b£\k) = e k ' 2 $_(k) + e 3k/ %(k) - e k ' 2 ^ + (k) - e k ^ 2 x[^ (k) , (75d 

b {2) (k) = e fe 0_ (k) - e- k l 2 xf ] (k) , (75e 

b (2 \k) = e k ^(k) + e%{k) - X^ik) - e~ k/2 X 2 [3) (k) , (75f 

b (2 \k) = e k ^(k) + e%(k) - qf\k) - AW(ife) , (75 g ; 

6f (jfe) = e fc ^_ (jfe) + e 2fc gi (k) - qf ] (k) - A® (jfe) , (751i 
bf \k) = e k 0_(k) + e 2 %(k) + e%(k) ~ e%(k) ~ x£\k) - e k ' 2 xf\k) , (75i 

6f(fe) = e fe ^(fc) + e 2k q 2 (k) ~ e k / 2 x[ 2 \k) , (75y 

b {3) (k) = e 3k/2 ^ (k) - x[ 3) (k) , (75k; 

b (3) (k) = e 3k / 2 ^(k) + e k /%(k) - X[ 2 \k) , (751 
bf\k) = e 3k / 2 ^(k) + e 3k /%(k) + e- k ' 2 x£\k) 

_ e 3fc / 2 gi(fc) - e- fe / 2 gf } (£;) - K 2 \k) , ^ 

b {3) (k) = e 3fc/2 ^_ (jfe) + e 5k/ % (k) - e 5fc/2 0_ (jfe) - X {3) (k) . (75n 

For the uppercase functions (x) we find 

Sf) (jfe) = e~ fe / 2 A« {k ) - e- k / 2 X^ (k) , (76a; 

(jfe) = e ~ k / 2 X^ (k) + X {2) (k) - e- fc / 2 ^ } (jfe) - A^ 2 ) (jfe) , (76b; 

B^\k) = e- k ' 2 qf\k) + e k / 2 x[ 1 \k) - e k l%(k) ~ xf\k) , (76c 

B ( l\k) = e k/ %(k) + e k ' 2 A^(k) - e k/2 + (k) - e k/2 x{ 1} (k) , (76d 

B {2) (k) = e- k / 2 A^ (k) - e- k ' 2 X {2) (k) , (76e 

B {2) (k) = qi h) (k) + e- fe / 2 Xf ) (k) - (k) - e~ k / 2 X® (k) , (76i 

B {2) (k) = X™ (k) + X™ (k) - 4 h) (k) - A« (k) , (76 g ; 

B {2) (k) = e k / 2 x[ 3) (k) + e" fc / 2 X {3) (k) - §f } (k) - A^ (k) , (76h; 

B (2) (k) = qf ) (k)+e k/2 X {2) (k) - e%(k) - X^ (k) - e k l 2 X {3) (k) , (76i 

B (2) (k) = e%(k) + e k/2 A (2 \k) - e k ' 2 X (2) (k) , (76j 

B (3) (k) = e~ k ' 2 A^ (k) - X (3) (k) , (76k 

B (3) (k) = e- k / 2 q^ (k) + X[ 3) (k) - X{ 2) (k) , (76f 
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(fc) = X {2) (k) + Xf ] (k) - e 3k/ % (k) - e- k/2 qi h) (k) - A< 2 > (k) , (76m) 
B^\k) = e 3k/ %(k) + e k/2 A {3) (k) - e 5k / 2 ^(k) - xf \k) . (76n) 

The crucial point is the observation that the latter forms a system of 14 linear equa- 
tions, exactly as many as there are unknown functions besides Bj n \k). Therefore, 
we can solve CED to get %(k), q 2 {k), q 3 (k), A«(fc), A^ik), A^(k), q£\k), q^\k), 
X[ 1] (k), X^\k), X[ 2 \k), X [2 \k), xf\k) and xf\k) in terms of the auxiliary func- 
tions B^ n \k). Eventually, we substitute this result into (|75p and are left with a set of 14 

equations, in which only b^\k) and B^\k) appear. 

We combine the results from the cases k < and k > to get a system of equations 
valid for all i£R. We find the equations to be 



b(fc) = — iiV smh(k(3/N) V(k) + K(fc) • B(Jfe) , (77) 



where 



b = 






B = 






V = 


(v {1) 


yW y(2) 



«2) t(3) W) 

>e ) °i j • • • i w 4 



f ,5? ),..., B^ 



J 



y{2) ^(3) T 7(3) 

[4] ' K [4] ' • • • ' V [4] 



(78) 

(79) 
(80) 



The matrix K(fc) is exactly the Fourier transform of the kernel matrix (|3ip . and the 
function 

v U) (k) _ sinh([g-j]fc/2) 

is the Fourier transform of (|38j) . Only the first term on the right-hand side of (I77p 
contains the Trotter number N explicitly. To analytically perform the global Trotter 
limit N — > oo, we therefore just have to consider 

lim Nsmh(kp/N) = k/3 . (82) 

Next, we apply the inverse Fourier transform and an integration over x to equation (|77|) to 
eventually obtain the NLIE (|26p - (|38p . The missing integration constants are determined 
by considering (|26|) in the limit x — > do. The asymptotics of the auxiliary functions can 
be easily read off from (|73p and (|74l) . because only the constant factors on the right-hand 
sides survive for large x. For the convolutions with the kernel functions, we find 

lim [Ki * B (n) ] (x) = £ (n) (oc) [°° Ki{x) ^ = « (n) (oo)^(0) • (83) 

x^oo L J 3 y_ oc 27T J 

Inserting this information into (|26[) finally leads to the constants given in (I39p . 
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To also derive the formula for the largest eigenvalue (|4U|) . one has to recall that the 
eigenvalue already appeared during the calculation. Its logarithmic Fourier transform 
A^(k) could be expressed solely in terms of Bj n \k). We additionally define the function 

_ (/>-(x - i)<f> + (x + 1) 

which has the advantage of having constant asymptotics. In the Trotter limit N — > oo, 
this definition simply yields lnA^)(0) = lnA^^O) — (3. Using the previously obtained 
information on A^(k), we can write 

A (1 \k) = iiVsirm(WiV)eHW2 g^Wg) + • lnB< B >(Jfc) . (85) 



Now we proceed as above. To perform the Trotter limit, we just have to use ([82]). Then 
we apply the inverse Fourier transform and an integration over the spectral parameter. 
Again, we find the integration constant by considering the limit x —* oo. We finally 
arrive at equation (|40p . 
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